Local stress and superfluid properties of solid 4 He 
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More than half a century ago Penrose asked [lt|: are the superfluid and solid state of matter 
mutually exclusive or do there exist "supersolid" materials where the atoms form a regular lattice 
and simultaneously flow without friction? Recent experiments provide evidence that supersolid 
behavior indeed exists in 4 He - the most quantum material known in Nature. In this paper we 
show that large local strain in the vicinity of crystalline defects is the origin of supersolidity in 4 He. 
Although ideal crystals of 4 He are not supersolid, the gap for vacancy creation closes when applying 
a moderate stress. While a homogeneous system simply becomes unstable at this point, the stressed 
core of crystalline defects (dislocations and grain boundaries) undergoes a radical transformation 
and can become superfluid. 

PACS numbers: 



The first microscopic mechanism of supersolidity 
was proposed by Andreev, Lifshitz Q and Chester 
Q arguing that a dilute gas of vacancies could 
lower the energy of an ideal helium crystal by quan- 
tum mechanically delocalizing vacancies. At low 
enough temperature the vacancies then undergo 
Bose-Einstein condensation and give rise to a su- 
persolid - and indeed vacancies are required for 
supersolidity 0]. Although the detection of fric- 
tionlcss mass flow was anticipated to be straight- 
forward, all experimental attempts failed for over 
thirty years [B| until Kim and Chan reported the 
observation of an unexpected drop of the torsional 
oscillator period at low temperatures 0, 0| ■ Since 
the period is proportional to the square root of 
the moment of inertia, the origin of the period 
shift was attributed to the supersolid decoupling 
of the 4 He solid mass from the rotating walls. Sev- 
eral groups have now confirmed the basic effect 
[SB EI EH but it is still unclear what physics lies 
behind these observations, both theoretically and 
experimentally [l^ ]. 

First-principle calculations [l3, 14 1 rule out the 
possibility of supersolidity in ideal equilibrium 
crystals of 4 He. On the experimental side, no 
steady flow of mass throu gh t he solid phase was 
detected in several setups jlol llij . It is believed 
that the flow detected at the melting curve in Ref. 
17| was along liquid channels which form on the 
cell walls in contact with the grain boundaries, i.e. 
it was not related to the supersolid phenomenon 

The situation has changed radically with the re- 



cent observation of the liquid pressure equilibra- 
tion by flow through a solid in the "UMass sand- 
wich" experiment reported by Ray and Hallock flij 
away from the melting curve. The dependence of 
the flow on the pressure difference was character- 
istic of critical superflow and inconsistent with the 
flow of a viscous liquid. In this work, we explain 
that the origin of superfluidity in solid 4 He is large 
local strain in the vicinity of crystalline defects. 
Superfluidity along dislocations or grain bound- 
aries can quantitatively account for the superflow 
through 4 He crystals observed by Ray and Hallock 

In the past several authors have speculated 
about the possibility of the stress-induced super- 
solidity [3, [l^, [2l[, especially under hydrostatic 
decompression since quantum effects and vacancy 
derealization are expected to increase at lower 
densities. The phenomenology missed the attrac- 
tion between vacancies which destabilizes a dilute 
homogeneous gas of vacancies. Central questions 
were left unanswered : what type and strength of 
stress is needed to close the insulating gap in 4 He 
and is this stress realistic? 

Our approach to this inherently strongly cor- 
related problem is numerical. Feynman's path- 
integral formulation of quantum mechanics allows 
an exact mapping between the 4 He system and a 
system of world lines in four-dimensional space. 
The fourth dimension represents evolution in imag- 
inary time and its extent is h/ksT where h is the 
the reduced Planck's constant, fcg is the Boltz- 
mann constant and T is temperature. This sys- 



60 
50 
40 
30 
20 
10 

-10 
-20 



n c n f 


n m 

I 































0.025 



0.0287 0.0315 

-3\ 



n(A~' 



FIG. 1: (Color online) Extrapolation of vacancy (Ay) 
and interstitial (A/) gaps shows that the density cor- 
responding to closing the gap lies in the liquid phase. 
The data points as a function of density for the ideal 
hep solid are taken from Ref. |23| . By nvil^i) we de- 
note the threshold chemical potentials for doping the 
system with vacancies(interstitials). 



tem of world lines is simulated efficientl y b y the 
quantum Monte Carlo "worm algorithm" [22| . En- 
ergy gaps for vacancy and interstitial excitations 
are readily obtained from the exponential decay in 
imaginary time of the one-body Green's function 
of the system [23| ■ 

Hydrostatic decompression After initial 
attempts to detect the supersolid state in 4 He 
failed, it was suggested that a metastable super- 
solid can form in crystals decompressed below their 
melting density of n m = 0.0287 A -3 (the freezing 
density is nf = 0.0261 A -3 ). The idea turned out 
to be impossible to realize experimentally, and now 
we understand why it was implausible in the first 
place. In Fig. Q] we show the density dependence 
of the vacancy and interstitial gaps. Though the 
simulation data are available for densities n > n m 
they can be readily and reliably extrapolated to 
lower densities n < n m using the near perfect lin- 
ear density dependence. One can easily see that 
metastable hep crystals remain insulating all the 
way to liquid densities and even beyond. It it un- 
likely that the solid structure will survive for long 
at liquid densities, and neither will the possible 
supersolid phase at density n c ss 0.025 A -3 . The 
hydrostatic strain required to reach this density is 
about [n m /n c — l) 1 / 3 ps 13.5 %. Assuming solid 



compressibility at melting 2_4j the required under- 
pressure is close to -25 bar. 

The stability of the metastable supersolid phase 
depends on the sign of the effective interactions 
between vacancies since an attractive Bose gas is 
unstable: for the vacancy gas this means a col- 
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FIG. 2: (Color online) Vacancy dispersion relation of 
the lowest band along the reciprocal lattice directions 
p z and p x for a system of 12 x 12 x 12 particles cal- 
culated at a density n = 0.0292A -3 and temperature 
T = 0.2K. For the hep lattice with two atoms in the 
unit cell there are two hopping amplitudes, one in the 
basal plane (t±), and one (t 2 ), along the TA direc- 
tion. At low temperatures the Monte Carlo method 
will project out the lowest branch of the dispersion re- 
lation. 



lapse to the lower density liquid. To determine 
the sign and strength of the vacancy-vacancy in- 
teraction one has to know their mass to* and 
their pair correlation function v(r). In the ideal 
gas v(r) is enhanced at short distances by a fac- 
tor of two relative to the large distance limit. 
Correspondingly, for repulsive/ attractive interac- 
tions v(r) is suppressed/enhanced relative to the 
ideal gas behavior. If attractive interactions are 
so strong that vacancies actually form a bound 
state, then v(r) is enhanced exponentially, and 
the decay of the correlation function with distance 
can be used to determine the binding energy from 
v (r) tx exp(— 2\/m*E(,r). 

In Fig. [2] we present the vacancy dispersion rela- 
tion which is analyzed within the tight-binding ap- 
proximation of Ref. pHj , obtaining tunneling am- 
plitudes t z = 0.45(5) K and t± = 0.50(5) K and 
effective masses m* = 0.45(5) and = 0.42(5) 
in units of the bare 4 He mass. We find that the 
effective mass is about two times lighter than that 
of 4 He atoms, and nearly isotropic. There is a 
small difference with the variational calculation of 
Ref. [25| which reported to^ = 0.31 in the basal 
plane and to* = 0.38 along the TA direction. 

In Fig. [3] we show the vacancy- vacancy correla- 
tion function, decreasing exponentially with dis- 
tance. This is concrete evidence showing that 
interactions between vacancies are attractive and 
strong enough to form a bound state. From the 




FIG. 3: (Color online) The probability of finding two 
vacancies at a distance r. The initial decay is fit with 
the exponential function v{r) oc exp(— 2^J m* Ebr) with 
Eh = 1.4(5)K. The simulation was performed for the 
hep solid with 6x6x6 lattice points and 214 atoms at 
a density n = 0.0292A" 3 and temperature T = 0.5K. 
At large distances the probability reaches a plateau due 
to finite temperature effects. 

exponential fit we estimate the binding energy to 
be Eb = 1-4(5) K. We conclude that the hep struc- 
ture at n c is kinetically unstable against collapse 
to the liquid state, and a supersolid in a decom- 
pressed crystals is not possible in the homegeneous 
setup. As a side remark we note that the equi- 
librium density of vacancies at finite temperature 
n v = (TO*T/27r) 3 / 2 e- Al '/ T « 0.015 n e" Av / T in- 
volves a very small pre-cxponcntial factor in ad- 
dition to the exponential suppression, which rules 
out their possible role in thermodynamics at tem- 
peratures below IK. 

Anisotropic stress We now turn to the study 
of crystals with non-zero diagonal components of 
the strain tensor u xx , u yy and u zz = —u xx ~ u yy . 
The phenomenology of the strain induced super- 
solidity in the hep structure [2(| is based on the 
minimal energy density allowed by symmetry: e = 
(a z — a±)u zz \ip\ 2 for the anisotropic stress and 
e = (a z u zz + a±(u xx + u yy ))\ip\ 2 for hydrostatic 
stress, where ip is the superfluid order param- 
eter. In the simulations we find ( see Fig. 2]) 
that a strain u zz = 2u xx = 2u yy of about 10- 
12% is necessary to close the gap for vacancy for- 
mation, and at higher strain the hep structure 
collapses. Such a strain corresponds to a stress 

C 2Z = C 33 U ZZ + Cl 3 (U XX + Uyy) = (C 33 - Cl 3 )U ZZ 

of approximately 50 bar, hardly achievable under 
realistic experimental conditions. Our data shows 
that within error bars there is initially no depen- 
dence of the gap on anisotropic compression. Thus 
the linear coupling to the anisotropic strain is close 



FIG. 4: (Color online) Gaps for vacancy formation un- 
der anisotropic compression. An ideal hep solid of size 
6 x 8 x 8 at density n = 0.0292 A~ 3 is subjected to 
anisotropic compression along the z direction. The vol- 
ume is kept constant by dilatation in the basal plane. 
The temperature is T = 0.25K. 



to zero, a z ~ a±, and one has to go beyond linear 
theory to account for the observed effects, and the 
closing of the gap. 

Shear stress By symmetry there is no linear 
coupling between \ip\ 2 and shear strain character- 
ized by the off-diagonal components u ZXl u yx , u zy . 
The anticipated strain dependence of the gap is 
quadratic, Ay = A^'[l — (u 2 zx + u 2 y/u 2 ]. Our re- 
sults in Fig.[5]allow us to estimate the critical value 
of shear strain as u c sa 0.15. Using measured val- 
ues of the elastic modulus C44 f=a 120-130 bar [26| 
the corresponding critical shear stress, 2Ci4U c , for 
closing an insulating gap in the hep He is about 
35 bar. 

Superfluidity along crystalline defects 

Even though the homogeneous strain-induced su- 
persolid phase is unstable, it can form locally if 
non-uniform strain accumulating close to struc- 
tural defects exceeds the critical value and desta- 
bilizes parts of the crystal. The most probable 
candidates for such highly strained superfluid de- 
fects are dislocations jia. ]20ll2lT|. with edge dislo- 
cations attracting most attention in the past be- 
cause they produce strain coupling linearly to the 
superfluid order parameter. Contrary to expec- 
tations, the first evidence for superfluidity in the 
dislocation core was reported for screw dislocations 
oriented along the z-direction j27|. These defects 
are characterized by the non-zero values of u zx 
which can be estimated by dividing the modulus 
of the Burgers vector b z — ^/8/3a by twice the cir- 
cumference of the circle going through the atoms 
closest to the core, Ana/y/3. The estimated strain 
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FIG. 5: (Color online) Gaps for vacancy formation un- 
der shear stress. An ideal hep solid of size L x x 12 x 12 
at density n = 0.0292A" 3 is subjected to shear stress 
along the z direction by uniformly deforming the sam- 
ple so that the atoms at boundary x — aL x (with 
a = 3.645A) are shifted exactly by one lattice period 
c = \/8/3a along the c-axis in order to match the pe- 
riodic boundary con ditio ns. The strain introduced is 
defined by u zx — y/2/3/L x and, for the chosen sizes 
L x — 12 — 5, the strain ranges from 0.07 down to 0.15. 
The temperature is T = 0.2K. 



most quantum solid in Nature. The actual struc- 
ture of defects is essentially unexplored territory, 
not less interesting than the solid matrix they re- 
side in. Similar exciting phenomena occur at de- 
fects and interfaces in electronic systems [H, . 
In this work we determined the critical values of 
the strain which are required to destabilize the hep 
structure of 4 He by closing its insulating gap, and 
find that these thresholds arc small enough to be 
exceeded at the dislocation cores and grain bound- 
aries. Local stresses thus hold the key to superso- 
lidity of 4 He. 

This work was supported by the National Sci- 
ence Foundation under Grants Nos. PHY-0653183 
and PHY-065135, CDRF grant No. 2853, the Nat- 
ural Science and Engineering Research Council of 
Canada under research grant G121210893, and the 
Swiss National Science Foundation. Simulations 
were performed on Hrcidar (ETH Zurich) , Typhon 
and Athena (CSI), and Masha (UMass) Beowulf 
clusters. 



^ju\ z + Uy Z = l/(v / 27r) ?s 0.22 exceeds the thresh- 
old value of u c = 0.15 found above. This explains 
why the superfluid density of the screw dislocation 
involves nearly all atoms closest to the nucleus. 

An alternative mechanism to release large strain 
at the core is to split the core into two parts, each 
characterized by half the Burgers vector. Whether 
this actually happens depends on the dislocation 
type and the energy of the domain wall stretching 
between the cores (the energy of stacking faults in 
hep lattices is relatively low). For screw disloca- 
tions the core splitting mechanism is not effective. 

Similar considerations apply to edge disloca- 
tions and grain boundaries, where simulations have 
detected superfluidity in some randomly oriented 
grain boundaries [28(, but the majority of small 
angle boundaries and stacking faults were found 
to be insulating. Simple estimates of strain in the 
vicinity of the defect cores range from 0.05 to 0.2, 
depending on the geometry and an effective Burg- 
ers vector, which may be reduced if the split-core 
mechanism is effective. 

Conclusions There is growing experimental 
evidence for unexpected properties of solid He at 
low temperature, including superfluid mass flow 
through the solid. The observed dependence on 
sample history, growth conditions, annealing and 
cooling procedures, indicate that crystalline de- 
fects are important for our understanding of the 
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